Quantum Phase Transitions of Light in the Dicke-Bose-Hubbard model 
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We extend the idea of quantum phase transitions of light in atom-photon system with Dicke- 
Bose-Hubbard model for arbitrary number of two-level atoms. The formulations of eigenenergies, 
effective Rabi frequencies, and critical chemical potentials for two atoms are derived. With a self- 
consistent method, we obtain a complete phase diagram for two two- level atoms on resonance, which 
indicates the transition from Mott insulator to superfluidity and with a mean excitations diagram 
for confirmation. We illustrate the generality of the method by constructing the dressed-state basis 
for arbitrary number of two- level atoms. In addition, we show that the Mott insulator lobes in the 
phase diagrams will smash out with the increase of atom numbers. The results of this work provide 
a step for studying the effects with combinations of Dicke-like and Hubbard-like models to simulate 
strongly correlated electron systems using photons. 

PACS numbers: 42.50.-p, 05.70.Fh 



I. INTRODUCTION 

Quantum phase transition (QPT) is the phase tran- 
sition that can only be accessed at absolute zero tem- 
perature, by the change of an external parameter or a 
coupling constant driven by quantum fluctuations 
QPT has attracted intensive studies in interacting many- 
body problems, originally for strongly correlated elec- 
tronic systems in condensed matter physics [2[ , and more 
recently for weakly interacting ultra-cold atomic system 
0] . Typically, it is difficult to control and probe such ex- 
otic quantum phenomena in strongly correlated systems 
of electrons. Optical lattices, artificial man-made crys- 
tals by interfering laser beams, offer a versatile platform 
for studying trapped Bose gases. In this situation Bose- 
Hubbard model including on site two atoms interaction 
and hopping between adjacent sites is used for describing 
the many-body dynamics from a Mott insulating to a su- 
perfluid phase in a gas of ultracold atoms with periodic 
potentials [J]. 

Instead photons are non-interacting Bosons, and there 
is no possibility to have any quantum phase transitions 
in purely photon systems. For pure Bose system, the 
conducting phase at zero temperature is presumably al- 
ways superfiuid f5[. However, engineered composites 
of optical cavities, few- level atoms, and laser light can 
form a strongly interacting many body system to study 
the concepts and methods in condensed matter physics 
from viewpoint of quantum optics. In this case, pho- 
tonic condensed-matter analogue is possible realized with 
state-of-art photonic crystals embedded with high-Q de- 
fect cavities. Therefore, photons interacting with atoms 
should be much easier to study and probe the critical 
quantum phenomena such as QPTs in conventional con- 
densed matter systems [fj. The simplest light-atom sys- 
tem is photons interacting with a single two-level atom 
(TLA), described by the Jaynes-Cummings model @, Q. 
With an array of high-Q electromagnetic cavities each 
containing a single TLA in the photon-blockade regime, 



quantum phase transitions of photonic insulator (excita- 
tions localization) to superfiuid (excitations derealiza- 
tion) are predicted by the Bose-Hubbard model @, [To[ 
and the XY spin model jllj . 

As the number of TLAs increases, collective effects due 
to interactions of atoms among themselves give rise to in- 
triguing many-body phenomena Q. In quantum optics 
Dicke model describes the collective spontaneous emis- 
sion of an initially excited ensemble of N TLAs interact- 
ing with a common photon field [l2j], and has triggered 
numerous investigations of various physical effects. A col- 
lection of atoms prepared in a certain initial state could 
decay collectively like a huge dipole with the emission 
of radiation not proportionally to the atom number N 
but to TV 2 , with a phase transition between a normal to 
a superradiant state [13j. Actually, when the maximal 
distance between any two of TLAs is much less than a 
typical wavelength, the coupling interaction for the pho- 
ton field no longer depends on the individual coordinate 
of the atoms but on the collective pseudo-spin coordi- 
nate. One has to add these pseudo-spins of N TLAs up 
to a single large pseudo-spin which is described by the 
collective angular momentum operators [8[ . 

Dicke model itself can offer the possible insights in the 
nature of QPT. With a size-consistent Hamiltonian for 
the Dicke model, squeezing of the photon field carries 
signatures of the associated quantum critical phenom- 
ena [lj]. QPTs of the Dicke-like and Bose-Hubbard-like 
models are both well-studied in optical and condensed 
matter systems, but the combination of these two mod- 
els are never studied. In this work, we extend the idea 
of QPTs of light in atom-photon system, proposed by 
A. D. Greentree et al. [9j], from a single TLA interac- 
tion system described by the Jaynes-Cummings model 
to arbitrary number of TLAs by the Dicke model. The 
problem we address here is TV identical TLAs which cou- 
ple to a single mode quantized radiation field within ideal 
photon cavities in the photon blockade regimes. With a 
self-consistent method, we numerically demonstrate that 
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Mott-insulator to superfluid quantum phase transitions 
exist even in the Dicke model for arbitrary number of 
two-level atoms interacting with photons. Detailed cal- 
culations of system eigenenergies, effective Rabi frequen- 
cies, and the critical chemical potentials for two TLAs 
are derived. The results of this work provide a more gen- 
eral picture for simulating strongly correlated electron 
systems using photons. 

This work is organized as follows, in Section [III we de- 
scribe the Dicke-Bose-Hubbard model used for studying 
QPT of light. The eigensystem solutions based on the 
dressed-state bases for two atoms are derived in Section 
Hill Results of mean field phase diagram, average exci- 
tations, and the extension to arbitrary number of TLAs 
are given in Section [iVl Section [V] gives the conclusion. 

II. DICKE-BOSE-HUBBARD MODEL 

Followed by the proposal by Greentree et al. [9] , Bose- 
Hubbard model can be extended to realize the Mott in- 
sulator to superfluid phase transitions in 2D photonic 
bandgap cavity network by including atom-photon inter- 
action. The Hamiltonian for our extended Dicke-Bose- 
Hubbard model is given by combining photon hopping 
between identical cavities in the photon-blockade regimes 
and the repulsive on-site TLAs interaction with the pres- 
ence of one-site energy chemical potentials, 

H = HP M ~ a+aj - (1) 

i ij i 

where z,j are the index for the individual cavities and 
range over all sites, Ni is the total number of atomic 
and photonic excitations. The second and third terms in 
the Hamiltonian, Eq.([T]), are Bose-Hubbard-like Hamil- 
tonian. The conserved particles in our system is the on- 
site total excitations Ni = af di + JfJ^. This conserved 
quantity are not pure photons, but the dressed photons 
which are mixtures of atoms and photons. We have as- 
sumed that all the inter-cavity hopping energy of photons 
K,i = Ki , and the chemical potential in the grand canoni- 
cal ensemble \±i=\± have no difference between cavities. 
The implementation of such photonic condensed-matter 
analogue is original proposed in Ref. Q . The first term in 
Eq.([T]) is the onset Hamiltonian for N TLAs interacting 
with a single mode field within a photon cavity, given by 
the Dicke model, i.e. 

HP M = eJ+ J- + Loa+ ai + /3(a, J+ + a+ jr), (2) 

where e is the transition energy for the TLAs, uo is the 
radiation field frequency, af and ai are the photon cre- 
ation and annihilation operators, and jf — ^ af, J~ = 
Y^j are the collective raising and lowering angular mo- 
mentum operators, respectively. The cavity mediated 
atom-photon coupling energy f3 is assumed to be real 
here. A superfluid order parameter 0, to be physically 
real (i.e. -0* = ifi), has been introduced for the studying 



of the QPTs in our Dicke-Bose-Hubbard Hamiltonian. 
The system is in the superfluid phase for a non-zero or- 
der parameter ip ^ 0, while in the insulator phase for a 
zero order parameter -0 = 0. With the mean field as- 
sumption (af) = V?*, we take decoupling approximation 
to investigate our Dicke-Bose-Hubbard Hamiltonian, i.e. 
afaj = ip*dj + ipaf — 1 0| 2 which is proportional to ip. 
Then the on-site mean field Hamiltonian of Eq.([T]) reads 

H™ F = H? M -KiP(ai + af)+ K \iP\ 2 -[iNi. (3) 

This mean field Hamiltonian is the same for every site. 
Without bothering the number of nearest neighbors 
around each idea photon cavity, we select three near- 
est neighbors per cavity in our simulation for the rea- 
son that it dose not affect our numerical results actually. 
The interaction part (3(diJf + afj~) in Eq.® actually 
does not change the field energy and commutes with ar- 
bitrary functions of the photon-number operator. Hence, 
we choose the eigenstates of the total excitations Ni to be 
the bare-states for the Dicke-Bose-Hubbard Hamiltonian. 
Consequently the subspace of Ni excitations spanned by 
Ni + 1 vectors is given by the direct pro duct of atom and 
field states, i.e. \atom)\photon) [HI, [l^. 

In the following, we choose two TLAs as an exam- 
ple, N = 2, for a clear illustration. The extension 
of the Dicke-Bose-Hubbard Hamiltonian for arbitrary 
number of TLAs is given in the Section IIVI which can 
be easily calculated by the same approach. For two 
TLAs, the bare-states are |0, e m )\n - 2), |#,e)|n-l), 
and |#® 2 ,0)|n) with photon number n which runs from 
0, 1, 2, 3 to oo. In our notation the collective angular mo- 
mentum eigenstates that describe the two TLAs are de- 
noted as |0, e® 2 ) for the case all the two atoms are in the 
excited state, |g, e) for the case only one atom is in the 
excited state, and |g® 2 ,0) for the case all the two atoms 
are in the ground state. These three bare-states are the 
normalized symmetric eigenstates of the noninteracting 
part ejfj~ -\-uafai in Eq.([2|). We neglect the dipole- 
dipole interaction, jf J J ,z ^ j, and represent the two 
TLAs interacting with the same cavity field simultane- 
ously in an idea photon cavity. For arbitrary excitations, 
in order to calculate the transition amplitude of Eq.®, 
we use these three complete symmetry degenerate bare- 
states replete with photons as the bases, i.e. 

|0,e® 2 )|0),| 5 , e )|l),|^ 2 ,0)|2), 
|0,e® 2 )|l),| 5 , e )|2),| 5 ® 2 ,0)|3), 

%e® 2 )\k-2),\g,e)\k-l),\g® 2 ,0)\k), 

|0,e® 2 )|n-2),| 5 ,e)|n-l),| 5 ® 2 ,0)|n). 

Here totally 3n bare-state bases form a group for the 
whole Hilbert space. Based on these on-site bases, we 
construct a 3n x 3n transition amplitude matrix for 
Eq.®, i.e. 



3 



2e - 2/i y/2/3 
^2/3 2e + uj - 3/i \/4/3 
Vi/3 2oj - 2/i 












-Kip 






-y/3Kijj 



2e + a; - 3/i V4/3 

\/4/3 2e + 2a; - 4/i ^6/3 
ygg 3a; - 3/i 



-V2ki/j 




(4) 



Eq.(|4j) is the starting matrix elements for our Dicke- 
Bose-Hubbard Hamiltonian with mean field assumption. 

III. EIGENSYSTEM SOLUTIONS 

In general it is very difficult to diagonalize Eq.(|4j) and 
find out all the desired dressed- states even for two TLAs. 
However, it may be instructive to diagonalize part of 
the Hamiltonian by assuming = for the mean field 
Hamiltonian. And a self-consistent method is applied for 
the case of ip ^ numerically. In Eq.(j4]), those entries 
in the boxes correspond to the block diagonal form of 
Eq.© with ji = therein, and the numbers of excita- 
tions in each block increases by one. As a matter of fact, 
we use the entries of the last box to calculate the cor- 



responding eigenenergies and eigenstates. Nevertheless, 
the analytic eigensystem solutions for non-zero detuning 
are very complicated, and we now focus on the case of 
zero detuning. In the follows we assume the system is on 
resonance, i.e. the atomic frequency is the same as the 
field frequency e = uj. In an idea photon cavity, when 
the two TLAs are excited, the rest of the photons will 
be dressed by the two TLAs. The interaction, last term 
in Eq.([2]), couples the three bare-states in the same ex- 
citation n manifold. The eigenstates and eigenvalues are 
derived in Eq.([5]|8]) for the center (£?| , n ), |0,n)) and up- 
per/lower {E\±^ n }, |=L,n)) branches. The eigenspectrum 
splits naturally into three branches, corresponding to the 
upper branch ^|+, n ), centre branch i£|o, n )> an d the lower 
branch E\_^ n y 



UJ 



E\o,n) = n- 



|0,n) 



|±,n> 



|±,n> 



-V^T|Q,e^ 2 )|n - 2) + y^ 2 ,0)|n) 
y/2n-l 
(2n + l)|±it>(n,|) 

2 

Vn|0, O \n - 2) + ^ [%±R(n, f )] \g, e) \n - 1) + V^T\g® 2 , 0) \n) 
2n-l + {^[*±R(n,*W 



(5) 
(6) 

(7) 
(8) 



FigtU shows the eigenenergies spectrum for two TLAs 
with the extended Dicke-Bose-Hubbard Hamiltonian. 
The interaction part, ^(a^J z + + J~), in Eq.© leads 
the 3n bare-states to couple together and form the nor- 
malized dressed- states. These three dressed- states for the 
n = 1 excitation manifold in Eqs.([5][8|) can be reduced 
to the so called triplet states and the isolated singlet 
state, |0,0) (l7[. Our definition for the dressed-states 
extend to n — is shown in FigHJ and we define the 
ground state for the dressed-state system as |0,0) with 
^|o,o) = 0- With a non-trivial form of the raising opera- 
tor, the ground state for the Dicke Hamiltonian is qual- 
itatively different from other dressed-states in Eqs.(|6]). 



The branches emerge at large resonant frequency, and the 
splitting increases with larger photon number, n, given 
by the effective Rabi frequency. 

To have a clear picture of the photon-atoms interaction 
in our two TLAs system, we introduce an effective Rabi 
frequency as R(n, ^) in Eq.(08|), which for n photons 
has the form, 

i?(n,^) = ^8(2n-l) + (^)2. (9) 

A two-state system has two possible states, i.e. a pho- 
ton in our model can either be in the excited or ground 
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state. The effective Rabi frequency induced here can 
be used as a measurement for the energy splitting be- 
tween the two states. For different photon numbers, the 
dressed-states described in Eqs.dSE]) oscillate with a dif- 
ferent Rabi frequency which is not only proportional to 
the photon number of the field but also with a depen- 
dence on the dimensionless resonant frequency uo/f3 as 
well. The dependence of the effective Rabi frequency on 
the photon number n for different normalized resonant 
frequency co//3 is shown in Fig|2j In the central region on 
the diagram, the boundary between domains defined by 
the effective Rabi frequency for different photon number 
eigenstate are the first state | — ,1), higher order states 
| — , 2) ,| — , 3), and so on. 

With the dressed-state formalism in Eqs.(|6ll8|), we can 
explain the features of Mott insulator to superfluid QPTs 
by a simple explanation [18]. Under the action of the 
Hamiltonian and the atomic population inversion oper- 
ator, the whole Hilbert space splits into one- and two- 
dimensional subspaces with n photons which are decou- 
pled. The probability amplitude achieves a saturate value 
when the number of photon increases. These dressed 
states are time independent and with constant ampli- 
tude, therefore when the atom- field system is prepared 
in a dressed states the atom remains stationary [l9j]. 
Since we know these dressed-states and their correspond- 
ing eigenenergies in Eqs.([5]|8]), the corresponding wave- 
function and the dynamics of the system will be just 
the superposition of them in addition to a phase term 
exp(— z£ , |{ ,±},n>)^ I n order to determine the ground 
state, we assume that ^|-, n ) < ^|+,n)> i- e - the nega- 
tive branch has lower eigenenergy. The three degener- 
ate bare-states for the n excitation manifold form the 
three nondegenerate dressed-states at three energies, i.e. 
n% - \ [R(n, !)_!],„" and nf + | [R(n, f ) - f ] . The 



principal feature of the Rabi frequency spectrum in Fig. 
[2]is easily understood in terms of the partly dressed states 
|+,n) and \ — ,n) in EqJH] which are separated in energy 
by an amount HR(n, ^). On the other hand, the cen- 
ter branch appears to be a number of the two branch 
£| +5n _i),£j_ n ) asymptotically. 

To derive analytical solutions for the dimensionless 
critical chemical potential, we substrate the negative 
branch energies, i.e. i£i n+ i) — E\_^, for that the sys- 
tem will change from n to n + 1 excitation per site and 
have the formula for different photon numbers, i.e. 



( UJ \ 
Mc(n, -) 



(10) 

V*fl(n, f )] 



2^2n(n - 1) 

As one can see in Fig. [3j the critical chemical poten- 
tial saturates to a constant value as the normalized res- 
onant frequency and the photon number increase. The 
dependence of the critical chemical potential with photon 
number is shown in the insert of Fig. [3l which is sensi- 
tive for the case of the normalized resonant frequency 
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FIG. 1: Eigenspectrum for two TLAs in the Dicke-Bose- 
Hubbard model, as a function of the normalized resonant fre- 
quency u)/ (3 for different photon numbers. The eigenspectrum 
splits naturally into three branches, which are shown as pos- 
itive i£|+,n) (upper), center E\ ^ (centered), and negative 
E\_^ (lower) with the corresponding dressed-states, |+,n), 
|0,n), and | — , n) , respectively. 
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FIG. 2: Effective Rabi frequency R(n,v/f3) defined in Eq.© 
v.s. normalized resonant frequency cu/f3. 



around w//3 « 10. This should be a suitable parameter 
for experimental observation [2(j. In order to calculate 
the ground state wavefunction of our generalized Dicke- 
Bose-Hubbard Hamiltonian we numerically solve Eq.(|4j) 
by applying a self-consistent method. When the photon 
number of the system increases, the ground state energy 
Eg will converges to the true ground state of the system. 
In Fig. [H we show the typical convergence of the ground 
state energy by increasing the photon numbers, which 
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Resonant frequency, go/(3 



FIG. 3: The critical chemical potential fi c (n,uj/ f3) v.s. nor- 
malized resonant frequency u>//3. The value of \i c saturates 
as oo /p and n increase in the limit with zero hopping k = 0. 
Insert shows the dependence of the critical chemical potential 
for different photon numbers. 



-10.5 




Photon number, n 



FIG. 4: The convergence of the ground state energy for dif- 
ferent photon numbers. 
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Photon hopping* *Yp 

FIG. 5: Phase diagram of the mean field Hamiltonian for the 
ground state on resonance, shown by the normalized inter- 
cavity hopping energy of photons k/(3 and the relative chem- 
ical potential (/i — co)//3. The notation SF refers to a super- 
fluid phase with strong interaction of photon hopping while 
MI refers to a Mott insulator phase with equally number of 
photons in each cavity. In the insulator region 1 0, 0) , | — , 1), 
and | — , 2) denote the negative branches of the dressed-states 
where the system will change from n to n + 1 excitation per 
site simultaneously filling photons in cavities and resulting in 
a finite gap of spectrum. 



2 excitation 




FIG. 6: The mean excitations for two TLAs. First few 
plateaus indicate constant density regions of the excitations 0, 
l,and 2 which correspond to the ground state configurations 
|0,0), | — ,1), and | — ,2), respectively. 



IV. RESULTS 



would quickly reach the minimum ground state energy 
for photon number up to n = 30. When stimulated with 
light, the interaction of each cavity with an ensemble 
of these atoms gives rise to a composite optical-atomic 
state as in Eqs.([5]l8j). However, if the photon number n 
is smaller than or just equal to the atom number TV, the 
quantum phase transition of light from Mott insulator to 
superfluid does not occur. Instead the ground state of the 
Dicke model exhibits an infinite sequence of instabilities 
quantum-phase-like transitions [15|. 



A. Mean field phase diagram 

After deriving the formulations for two TLAs in our 
system, we calculate the phase diagram for the Dicke- 
Bose-Hubbard Hamiltonian in Eq.(|3]) for two TLAs by 
applying the mean field theory and the self-consistent 
method, as shown in Fig|5l Clear quantum phase transi- 
tions can be seen in the diagram for different normalized 
inter-cavity hopping energy of photons k,//3 and different 
relative chemical potential (/i — uo)/ f3. Here the notation 
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(a) 3 TLAs 



(b) 4 TLAs 





FIG. 7: The phase diagrams for arbitrary number of TLAs, 
(a) N = 3, (b) N = 4, (c) N = 5, (d) N = 6, (e) N = 7, and 
(f) N = 10. The axises of horizontal and vertical are defined 
as the same as in Fig(5] The phase boundary between MI-SF 
are the superposition states that disappear and result in a 
single macroscopic coherent radiation state for large number 
of TLAS. 



SF refers to a superfluid phase with strong interaction of 
photon hopping. And the notation MI refers to a Mott 
insulator phase with equally number of photons in each 
cavity. The Mott insulator to superfluid phase transitions 
occur for the case of the photon number is larger than the 
atom number, N <C n. We give a simple picture for the 
QPTs of light in our system. As photons pass through an 
array of high-Q cavities with two TLAs per site, there is 
an upper limit for the energy of the two TLAS as all the 
atoms have been excited to the upper level. At the end, 
interaction photons are not independent to each other 
due to the iterations from effective on-site repulsion, and 
inter-cavity hopping with their proximity. In such a way 
the two TLAs are strongly correlated. Change of the 
critical chemical potential results in the change of the 
on-site numbers of photons, as shown in FigEl Regard- 
ing to the on-site chemical energy /i, the regions to the 
right corresponding to chemical potential \i ^ In 
the insulator region |0,0), |— ,1), and |— ,2) denote the 
negative branches of the dressed-states where the system 



will change from n to n + 1 excitation per site, simulta- 
neously filling photons in cavities and resulting in a finite 
gap of spectrum. The superfluid phase is the eigenstates 
of cii and excitations over the |— , n) branches in the FigJSJ 
The probability of finding the average photon number n 
in this regime obeys Poisson distribution. With increas- 
ing laser intensity n <C /3, we have strong interactions. It 
is the most important regime with rich dynamics where 
the on-site repulsion dominates with equally numbers of 
photons in each cavity as one can see on the left in FigEJ 
This region is corresponding to a constant density of pho- 
tons filling in cavities simultaneously. In this situation 
each site has exactly the same integer number of photons 
with strongly on-site coupling regime. A finite gap of 
spectrum is formed and photons here are incompressible, 
resulting an insulator phase on the other hand H, d, [2l| . 



B. Average excitations 

Another way to indicate QPTs in our Dicke-Bose- 
Hubbard Hamiltonian is to study the average excitations. 
We minimize the ground state energy with respect to 
the order parameter, -0, and take the derivative of the 
ground state energy with respect to p. Theoretically, we 
can consider the problem by fixing the chemical poten- 
tial by varying photon numbers in the grand-canonical 
ensemble and the derivatives value can be expressed by 
the ensemble average, i.e. 



dp 



\lp = 1prr 



(11) 



The mean excitations for two TLAs in our systems are 
shown in Fig. [6l It can be clearly seen that the density of 
the incompressible Mott insulator phases does not change 
with the change of the relative chemical potential, and 
the value of the relative chemical potential jumps discon- 
tinuously goes through a lobe. Regions with varying p 
have coherent states as ground state configurations. In 
fact, we use the mean excitations to confirm the num- 
bers of photons in each Mott lobe, and it is the mean 
excitation p. On can view this average excitations for 
the evidences of the phase transitions and lobes manifest 
in Fig. [5l 



C. Extension to arbitrary number of TLAs 

For arbitrary number of TLAs, A/", and arbitrary num- 
ber of photons, n, we use following general bases for 
Eq.(j3]) to solve our Dicke-Bose-Hubbard Hamiltonian, 
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\o,e® N )\o),\g® {N - 1 \e}\i),---\ g ,e®( N -V)\N-i),\g® N ,o)\N}, 

|0, e ^)|l),|^- 1 ), e )|2>,...| 5 ,e^- 1 ))|iV),|^,0)|iV+l), 

|0, e® N ) \k-N), Ig^-'le) \k - N + 1), • • • \g, e®^" 1 )) \k + N-1), \g® N , 0)\k), 

|0, e 0iV ) |n - N), Ig^-V, e) \n - N + 1), • • • \g, e®^" 1 )) |n + N - 1), 0)|n). 



Strong coupling with many excitations is expected to be 
first observed for a small number of N [22[ , there we pro- 
vide numerical results for a comparison from smaller to 
larger value of atom numbers N. By taking the limitation 
of an arbitrary numbers of TLAs with fixed photon num- 
bers up to 30, we obtain the Mott insulator to superfluid 
phase diagrams for the number of TLAs from N = 3 to 
N = 10 in FigJTl respectively. By increasing the number 
of TLAs with the same amount of deposited photons, 
resulting in a insulating to superfluid phase transition 
which is characterized by the MI-SF transitions of regu- 
lar bosons on the regions above the first lobe. The value 
of k at the tip of the n-th Mott lobe varies as ~ for 
large N. There is no energy barrier to the addition of ex- 
traphotons and superfluidity occurs at arbitrarily small 
k, [5|. When the limits N — » oo or (3 — » oo or both are 
taken, one has a single macroscopic coherent radiation 
state and the superposition of dressed-states disappears 
[23( ] . We can thus conclude that as the number of TLAs 
increases, superposition states may be disappearing and 
classically emerges. 



transitions with photons can be realized in an extended 
Dicke model for arbitrary number of two- level atoms. We 
illustrate the generality of the method by constructing 
the dressed-state basis for arbitrary number of two-level 
atoms. Moreover, we show that as the number of TLAs 
increases, superposition states may be disappearing and 
classically emerges. With more controllable lightwave 
technologies, the understanding of quantum phases tran- 
sitions of light with distinctive properties, organizations 
of the ground state wave function, and practicable new 
applications should make it more easily to be realized. 
With combinations of Dicke-like and Hubbard-like mod- 
els to simulate strongly correlated electron systems us- 
ing photons, we believe that there would be more and 
more interesting quantum phase transitions of light to 
be demonstrated as those in condensed matter physics. 



V. CONCLUSION 

With the Dicke-Bose-Hubbard Hamiltonian, we show 
that the Mott insulator to superfluid quantum phase 
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